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Abstract. We consider the random matrix ensemble with an external source 

defined onnxfi Hermitian matrices, where A is a diagonal matrix with only two eigenvalues 
±a of equal multiplicity. For the case a > 1, we establish the universal behavior of local 
eigenvalue correlations in the limit n — *■ oo, which is known from unitarily invariant random 
matrix models. Thus, local eigenvalue correlations are expressed in terms of the sine kernel 
{SJ , in the bulk and in terms of the Airy kernel at the edge of the spectrum. We use a character- 

ization of the associated multiple Hermite polynomials by a 3 x 3-matrix Riemann-Hilbert 
problem, and the Deift/Zhou steepest descent method to analyze the Riemann-Hilbert prob- 
lem in the large n limit. 
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1. Introduction and statement of results 
We will consider the random matrix ensemble with an external source, 



defined onnxn Hermitian matrices M. The number n is a large parameter in the ensemble. 
The Gaussian ensemble, V(M) = \M 2 1 has been solved in the papers of Pastur [21] and 
Brezin-Hikami [7| ^Uj> by using spectral methods and a contour integration formula for the 
determinantal kernel. In the present work we will develop a completely different approach 
to the solution of the Gaussian ensemble with external source. Our approach is based on 
the Riemann-Hilbert problem and it is applicable, in principle, to a general V. 

We will assume that the external source A is a fixed diagonal matrix with n\ eigenvalues 
a and n 2 eigenvalues (—a), 

(1.2) 




As shown by P. Zinn- Justin j^Zj, for any m > 1, the m-point correlation function of eigen- 
values of M has the determinantal form, 

l<j,k<m- 

(1.3) 
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In our previous work [E] we show that the kernel K n (x,y) can be expressed in terms of a 
solution to the following matrix Riemann-Hilbert (RH) problem: find Y : C \ R — > C 3x3 
such that 

• Y is analytic on C \ 1, 

• for x G R, we have 

(1 wi(x) w 2 (x)\ 
1 , (1.4) 
1 J 

where 

Wl (x) = e - n(y{x) - ax \ w 2 (x) = e -"(n*)+™) y (1.5) 

and Y + (x) (Yl(x)) denotes the limit of Y(z) as z — > x from the upper (lower) half- 
plane, 

• as z — > oo, we have 

y w=(' +o G))(s v JLJ- (i - 6) 

where / denotes the 3x3 identity matrix. 
Namely, 

]e nay[ Y (y)-'Y( x)} 21 + e- na y[Y(y)- l Y(x)}^ 
2iri(x — y) 

(1.7) 



e -Jn(V(*)+V(y)) 
27TZ(X — ?/) 



(0 e nay e~ nay ) Y{y)' 1 Y{x 




The RH problem has a unique solution and the solution is expressed in terms of multiple 
orthogonal polynomials, see |H] and Section 2.1 below. For now, let us mention that the 
(1, 1) entry Y\\ satisfies 

Y u (z) =E[det(zI - M)} (1.8) 

where E denotes expectation with respect to the measure (jl.lj) . So it is the average charac- 
teristic polynomial for the random matrix ensemble. 

It is the aim of this paper to analyze the RH problem as n — ► oo, by using the method of 
steepest descent / stationary phase of Deift and Zhou jT3] . We focus here on the Gaussian 
case V(x) = \x 2 . Our first result concerns the limiting mean eigenvalue density. 

Theorem 1.1. Let V(M) = \M 2 , n\ = n 2 = n/2 (so n is even) and let a > 1. Then the 
limiting mean density of eigenvalues 

p{x) = lim — K n (x, x) (1.9) 

n-^oo n 

exists, and it is supported by two intervals, \—Z\^—z 2 \ and [z 2 ,Zi\. The density p(x) is 
expressed as 

p(x) = -\lm£(x)\, (1.10) 

7T 
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where £ = £(x) solve the cubic equation, 

£ 3 -x£ 2 - (a 2 - l)£ + xa 2 = (1.11) 

(Pastur's equation). The density p is real analytic on (— z%, — z 2 ) U {z 2 ,zi) and it vanishes 
like a square root at the edge points of its support, i.e., there exist constants pi,p 2 > such 
that 

p(x) = — \x — z,-| 1//2 (l + o(l)) as x — * Zj, x G (z 2 ,Zi), 



"31 

p(x) = — \x + Z,-| 1/,2 (l + o(l)) as 2 ^ — Zj, x G ( — ^i, — z 2 ). 



71 

P± 
7Y 



1.12) 



Remark: We obtain p from an analysis of the equation 

_ e - (a 2 - i)e 
e - « 2 ' 

The critical points of the mapping ()1.13|) satisfy 

1 



;i.i3) 



f 2 = - + a 2 ±-Vl + 8a 2 . (1.14) 

For a > 1, the four critical points are real, and they correspond to four real branch points 
±zi, ±z 2 with z\ > z 2 > 0. We denote the three inverses of ()1.13j) by £j(z), j = 1,2,3, 
where £i is chosen such that £i(z) ~ z as z —>■ oo. Then £i has an analytic continuation to 
C \ ([— zi, —z 2 \ U [z 2 , Zi\) and Im£i + (x) > for x G (— Zj, —22) U (z 2 , z\). Then the density p 
is 

p(x) = -lm£ 1+ (x), (1.15) 

7T 

see Section 3. 

The assumption a > 1 is essential for four real branch points and a limiting mean eigen- 
value density which is supported on two disjoint intervals. For < a < 1, two branch points 
are purely imaginary, and the limiting mean eigenvalue density is supported on one inter- 
val. The main theorem on the local eigenvalue correlations continues to hold, but its proof 
requires a different analysis of the RH problem. This will be done in part II. In part III we 
will discuss the case a = 1. 



Remark: The density p can also be characterized by a minimization problem for logarithmic 
potentials. Consider the following energy functional defined on pairs (pi, p 2 ) of measures: 



E(pi,p 2 ) =11 log- rdpi{x)dpi{x) +11 log- ;dp 2 (x)dp 2 (y) 

JJ \x-y\ JJ \x-y\ 

+ J J log ^ 1 ^ dp 1 (x)dp 2 (y) + J Qx 2 - ax^j dpi(x) + J Qx 2 + ax^j dp 2 (x). 

The problem is to minimize E(px, p 2 ) among all pairs (pi, p 2 ) of measures on K with J dp\ = 
J dp 2 = |. There is a unique minimizer, and for a > 1, it can be shown that p\ is supported 
on [z 2 , Zi), p 2 is supported on \—z\, —z 2 ] and p is the density of pi + p 2 . This minimal energy 
problem is similar to the minimal energy problem for Angelesco systems in the theory of 
multiple orthogonal polynomials, see pfl I17j. 
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It is possible to base the asymptotic analysis of the RH problem on the minimization 
problem, as done by Deift et al, see ^21 E3 El, for the unitarily invariant random matrix 
model. However, we will not pursue that here. 

Our main results concern the universality of local eigenvalue correlations in the large n 
limit. This was established for unitarily invariant random matrix models 

by Bleher and Its for a quartic polynomial V, and by Deift et al ^3] for general real 
analytic V. The universality may be expressed by the following limit 

i- 1 T ^ ( u v \ sin7r(u — v) 
hm — - — -K n x H — r,a;oH -, — r = -, — (1.17) 

np(x ) V np(x ) np(x ) J ir(u-v) 



n^oo 



which is valid for xq in the bulk of the spectrum, i.e., for xq such that the limiting mean 
eigenvalue density p(xo) is positive. The proof of ()1.17|) established Dyson's universality 
conjecture pHl for unitary ensembles. 

In our case, we use a rescaled version of the kernel K n 

K n (x,y) = e n ^- h ^K n (x,y) (1.18) 

for some function h. The rescaling (jl,18|) is allowed because it does not affect the correlation 
functions R m (ll.Hj) . which are expressed as determinants of the kernel. Note that the kernel 
K n of (jl.7|) is non-symmetric and there is no obvious a priori scaling for it. The function h 
in ()1.18J) has the following form on (— z\, —z 2 ) U (z 2 , Zi) 

h{x) = — -x 2 + Re Ai+(x), x 6 (— z±, — z 2 ) U (z 2 , z\) (1-19) 
with Ai + (a;) = £i + (s)ds, where £i is as in the first remark after Theorem ll.il 

Theorem 1.2. Let V(M) = \M 2 , n x = n 2 = n/2, and let a > 1. Let Zi,z 2 and p be as 
in Theorem 11.11 and let K n be as in (I1.18j) . Then for every xq G (— z%, —z 2 ) U (z 2 , Z\) and 
M,oGl, we have 

hm ^—:K n x + — rT , x + —— = ^ (1.20) 

np{x ) \ np(x ) np(x )J ir(u-v) 



n— >oo 



Our final result concerns the eigenvalue correlations near the edge points ±Zj. For uni- 
tarily invariant random matrix ensembles ()1.16|) the local correlations near edge points are 
expressed in the limit n —>■ oo in terms of the Airy kernel 

Ai(u)Ai'(v) - Ai'(u)Ai(v) ^ 
u — V 

provided that the limiting mean eigenvalue density vanishes like a square root, which is 
generically the case [T^. In our non-unitarily invariant random matrix model, the limiting 
mean eigenvalue density vanishes like a square root, ()1.12|) . and indeed we recover the kernel 
(jl.21j) in the limit n — > oo. 
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Theorem 1.3. We use the same notation as in Theorem M. 21 Let p\ and p 2 be the constants 
from (jl.l2|) . Then for every u,v G 1R we have 

1 T y ( u v \ Ai(u)Ai» - Ai'(u)Ai(v) 
km — K n z x + ^ , Zl + ^ = K J 1 ; W , 1-22 

- {piny/ 6 \ (piny/ 6 \p\nyl A ) u — v 



and 

km 1 K (z U z V ^ - Ai ( M ) Ai '( t ') ~ A[ '( U ) A K V ) q 23 x 

Similar limits hold near the edge points —z% and —z 2 . 

As said before, our results follow from an asymptotic analysis of the RH problem (jl.4jl - 
(jl.fijl . which involves 3x3 matrices. In the past, asymptotics for RH problems has mostly 
been restricted to 2 x 2-matrix valued RH problems, see e.g. jH EH Clj and references 
cited therein. The first asymptotic analysis of a 3 x 3 matrix RH problem appeared in [20.. 
in an approximation problem for the exponential function. In the present work we use some 
of the ideas of |2U] . 

As in [201 a main tool in the analysis is an appropriate three sheeted Riemann surface. To 
motivate the choice of the Riemann surface we describe in Section El the recurrence relations 
and differential equations that are satisfied by a matrix ty, which is an easy modification of 
Y, see (12.7)1 below. The Riemann surface is studied in Section El and we obtain from it the 
functions £j and \j, j = 1, 2, 3, that are necessary for the transformations of the RH problem. 
The first transformation Y i— > T normalizes the RH problem at infinity and at the same time 
introduces oscillating diagonal entries in the jump matrices on the cuts [—zi, —z 2 ] and [z 2 , z\], 
see Section 0] The second transformation T i— > S involves opening of lenses around the cuts, 
which results in a RH problem for S with rapidly decaying off-diagonal entries in the jump 
matrices on the upper and lower boundaries of the lenses, see Section The next step is the 
construction of a parametrix, an approximate solution to the RH problem. In Section El we 
ignore all jumps in the RH problem for S, except those on the cuts [—Zi, —z 2 ], [z 2l Zi]. This 
leads to a model RH problem, which we solve by lifting it to the Riemann surface via the 
functions This leads to the parametrix away from the edge points ±zi, ±z 2 . A separate 
construction is needed near the edge points. This is done in Section [7| where we build the 
local parametrices out of Airy functions. The final transformation S \— > R is done in Section 
[HJand it leads to a RH problem for R whose jump matrices are uniformly close to the identity 
matrix. Then we can use estimates on solution of RH problems, see jT^j, to conclude that 
R is close to the identiy matrix, with error estimates. Having that we can give the proofs of 
the theorems in Section 

Our approach proves simultaneously large n asymptotics of the (1, 1) entry of Y, which by 
(|1.8|) is equal to the average characteristic polynomial. This polynomial is called a multiple 
Hermite polynomial for the case of V(x) = |x 2 , see [H] and Section El below. Since its 
asymptotics may be of independent interest, we consider it briefly in Section El below. 
More information on multiple orthogonal polynomials and their asymptotics can be found 
in [T71 122 122] , see also the surveys [TJ 0] and the references cited therein. 
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2. Recurrence relations and differential equations 

In order to motivate the introduction of the Riemann surface associated with (jl.l3j) we 
discuss here the recurrence relations and differential equations that are satisfied by the 
solution of the Riemann-Hilbert problem (|1.4|) - (|1.6j) in case V(x) = \x 2 . It also reveals the 
integrable structure. We note however, that the results of this subsection are not essential 
for the rest of the paper. 

For the recurrence relations we need to separate the indices n\ and n 2 in the asymptotic 
behavior from the exponent n in the weight functions w±, w 2 of (jl.5|) . In this section 
we put 

Wl {x) = e - N ^ x2 ~ ax \ w 2 (x) = e -N(^ 2 +ax) ^ 

where N is fixed, and we let Y = Y ni>m be the solution of the Riemann-Hilbert problem 
(jl.4j) . f)1.6|) with V(x) = \x 2 and w±, w 2 given by (j2.1|) . Let P ni;Tl2 (x) = x n + ■ ■ ■ be a monic 
polynomial of degree n = n\ + n 2 such that for k = 1,2, 

P ni ,n 2 (3o)x 3 w k (x)dx = 0, 0<j<n k -l. (2.2) 

The polynomial P nitm (x) is unique and it is called a multiple Hermite polynomial, see |2^l2*o]. 
Denote for k = 1,2, 

/oo 
P num (x)x n «w k (x)dx ^ 0. (2.3) 
-oo 

The solution to the RH problem is 

(Pm,n 2 C(P ni ^ n2 Wi) C(P ni ,n2 W 2) \ 

l,n 2 W X ) CiC{P ni ^ hn2 W 2 ) , (2.4) 

c 2 C{P ni , n2 ^iWi) c 2 C{P ni)n2 _ x w 2 )) 

with the constants 

2ni 2m 

c i--7(i) > c 2 --7(2) ' y 2 - b > 

" / ni-l,n 2 ll n 1 ,n 2 -l 

and where Cf denotes the Cauchy transform of /, 

Cf(z) = ^- [ WLds. (2.6) 
2m J R s- z 

The recurrence relations and differential equations are nicer formulated in terms of the 
function 



P ni , m e-* Nz2 C{P ni , n2Wl )e- Naz C(P nun2 w 2 )e Naz 

y ni ,n 2 = | Pni-i.nae-3"* 3 C (P^^w^e^ C{P ni ^ n2 w 2 )e Naz 

K P ni ,n 2 -ie~* Nz2 C{P ni , n ^ lWl )e- Naz C{P ni)n2 . lW2 )e Na - 

'10 

-i 




(2.7) 



q 1 | Y nur 
c^ 1 

The function \I/ = ^f ni ,n 2 solves the following RH problem: 
• \I/ is analytic on C \ 1, 
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for x G R, we have 



as z — > oo, we have 



'1 1 1 N 
*_(.<■) = #_(.,■) |010 
,0 1. 



(2.f 



i + o[ - 

z 










c ~l^— nig— Afaz 










(2.9) 



Proposition 2.1. VKe aai>e £/ie recurrence relations, 



^m,n 2 +l(^) 

and the differential equation, 



ni 


n-z 


N 


N 











-2a 


711 


ri2 




N 


2a 












(2.10) 





1 


"1 

iV 

—a 


"2 

iV 




1 





a 



(2.11) 



The proof of Proposition 12.11 is given in the Appendix [0 below. 

We look for a WKB solution of the differential equation ()2.11)1 of the form 

^ ni ,n 2 (z)=W(z)e- Nk ^\ 



(2.12) 



where A is a diagonal matrix. By substituting this form into ()2.11|) we obtain the equation, 



-WA'W' 1 = A WW' 1 , 

N 



(2.13) 



where A is the matrix of coefficients in ()2.11j) . By dropping the last term we reduce it to the 
eigenvalue problem, 

WA'W' 1 = -A. (2.14) 

The characteristic polynomial is 



det [£I + A] 



£ — z h t 2 

- 1 £- a 

- 1 i + a 

£ 3 - z£ 2 + (ti +t 2 - a 2 )£ + (ii -t 2 + za)a, 



(2.15) 



where t\ = ^ and t 



— B2. 

2 N ■ 



The spectral curve £ 3 — z^ 2 + {t\+t 2 — a 2 )£+ (t\ —t 2 + za)a = defines a Riemann surface, 
which in the case of interest in this paper (where N = n and n\ = n 2 = |n) reduces to 



^3 _ ^2 _ (fl 2 _ + za 2 



0. 



(2.16) 



This defines the Riemann surface that will play a central role in the rest of the paper. 



(3.i; 
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3. RlEMANN SURFACE 

The Riemann surface is given by the equation ()2.16|) or, if we solve for z, 

e - (a 2 - l)g 

e - ^ ■ 

There are three inverse functions to ()3.1|) . which we choose such that as z - 

«*) = »-i + o(l 

&(2)=o+ J_ + (^, (3.2) 

6W =-a + ± + (l 

We need to find the sheet structure of the Riemann surface (j2.16jh The critical points of 
z(£) satisfy the equation 

£ 4 -(l + 2a 2 )£ 2 + (a 2 -l)a 2 = 0, (3.3) 
which is biquadratic. The roots are 

£i 2 )2 = ^ + a 2 ±~vT+8^. (3.4) 

The value a = 1 is critical, in the sense that for a > 1 all the roots are real, while for a < 1, 
two are real and two are purely imaginary. In this paper we will consider the case a > 1. As 
noted before, we will consider the cases a < 1 and a = 1 in parts II and III. 
Set 

p,g=^ + a 2 =F^v / i + 8a 2 , < p < q. (3.5) 

Then the critical points are £ = ±p, ±q. The branch points on the z-plane are ±Zi and 
±2 2 , where 



VI + 8a 2 + 3 VI + 8a 2 - 3 

Zl = g— F= , 2 2 = P~F= j < z 2 < Zl- (3.6) 

VI + 8a 2 + 1 Vl + 8a 2 -l 

We can show that £i, £ 2 , an d £3 have analytic extensions to C \ ([—21,-22] U [22,21]), 
\ [2 2 , Zi] and C \ [—21, — 2 2 ], respectively. Also on the cut [z 2 , 21], 



£i+(z) = 6-(z) = 6-(^) = 6+ fa), ^2 < x < zi 
Im£i + fa) > 0, z 2 < £ < 21 

and £3(2;) is real. On the cut [—21, — z 2 ], 



(3.7) 



(3i 



6+ (a?) = 6- fa) = £3- fa) = 6+ fa), - zi < x < -2 2 , 
Im£i + (x) > 0, — zi < x < — z 2 , 

and £ 2 fa) is real. Figure 1 depicts the three sheets of the Riemann surface (I2.16JI . 
We define 

p(x) = -Im£i + (x), x e [-21, -2 2 ] U [z 2 , zj. (3.9) 

7T 
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Figure 1. The three sheets of the Riemann surface 



Proposition 3.1. We have p{x) > for x G (—Z\, —z 2 ) U (z 2 , z{) and 

r-z 2 rz! 1 

/ p(x)dx = / p{x)dx = -. 
J -zi J z 2 2 

Moreover, there are p\,pi > such that 



(3.10) 



p(x) = — \x — Zj\ 1/2 (1 + 0(x — Zj)) as x — > zj, x G (z 2 , z\) 

7T 

p(x) = — \x + Zj\ l l 2 (1 + 0(x + z,-)) asx-> — Zj, x G (— z\, —z 2 ). 

7T 



(3.11) 



Proof: The fact that p(x) > for x G (— 2i, — £2) U (-22, 21) was already noted in (|3.7|1 and 

(EH). 

We have for x G [^ 2 , ^1], 
Thus 

y = ^-j> ii{z)dz 

where 7 is a contour encircling the interval [z 2 , z±] once in the positive direction. Letting the 
contour go to infinity and using the asymptotic behavior (|3.2|) of £2(2) as z — > 00, we find 
the value of the second integral in f)3. lOj) . The value of the first integral follows in the same 
way. 

For (j3.11|) we note that near the branch point z\, we have for a constant p\ > 0, 



Uz) = q + Pi(z-z 1 ) 1 / 2 + 0(z-z 1 ) 
&(z) = q-pi(z-z 1 ) 1 / 2 + 0(z-z 1 ) 



(3.12) 
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as z — > z\. Similarly, near Z2 we have for a constant p 2 > 0, 

t-iiz) =p- p 2 (z 2 - zf' 2 + 0(z 2 - z) 

U*) = P + P2(?2 - Z) 1 ' 2 + O (Z 2 - Z) 

as z — > z 2 (with main branches of the square root). By symmetry, we have similar expressions 
near —z\ and — z 2 and (j3.11j) follows. □ 

Next, we need the integrals of the ^-functions, 

h(z) = f Us)ds, k = 1,2,3, (3.14) 



which we take so that Ai and A2 are analytic on C \ (—00, z±] and A3 is analytic on C \ 
(—00, —z 2 \. From ()3.2|) it follows that, as z — > 00, 

z 2 ( 1 

\ x {z) = — -Inz + h + O f — 

A 2 (z) = az + ^ln,2 + Z 2 + (- ) , (3.15) 



2 

1 /l 

A 3 (z) = -az + - Inz + k + O [- 

where /1, / 2 , £3 are some constants, which we choose as follows. We choose l\ and l 2 such 
that 

Ai(z0 = A 2 (^) = 0, 

and then Z3 such that 

A 3 (-^2) = Ai+(-2 2 ) = Ai_(-2; 2 ) - Txi. 
Then we have the following jump relations: 

A i+ (x) - Ai_(x) = —717, X G [—2 2 , 32], 

Ai + (x) — Ai_(x) = —2ni, x G (—00, —zi], 

A 2 +(s) - A 2 _(x) = 7rz, x G (-00, z 2 ], 

Ai+(x) = A 2 _(x), Ai_(a;) = A 2+ (x), xe[z 2 ,zx), 

A i+ (x) = A 3 _(x), Ai_(x) - 7rz = A 3+ (x), x G [-zi, -z 2 ]- 

A 3+ (x) — A 3 _(x) = 7T2, X G ( — OO, — 2],]. 



(3.16) 



Note that due to the first two equations of (J3.16|) we have that e nXl ^ is analytic on the 
complex plane with cuts on [—z\, —z 2 ] and [z 2 ,z\] (recall that n is even). Furthermore, we 
also see that e nA2 ^ (resp., e nXa ^) is analytic on the complex plane with a cut on [z 2 ,Zi] 
(resp., [—Zi,—z 2 ]), see Figured] 

For later use, we state the following two propositions. 

Proposition 3.2. On M \ [z 2 ,z\] we have ReA 2+ < ReAi_ ; and onR\ [— z±,— z 2 \, we have 
Re A3 + < Re Ai_. 
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Proof. It is easy to see that > ^(x) for x > z\. Since \\{z\) = \2{zi) and A^ = £j for 

j = 1, 2, 3, it is then clear that Xi(x) > \ 2 (x) for x > Z\. 

We also have that Re£i_(x) < £, 2 {x) for x < z 2 , from which it follows that ReAi_(x) > 
Re A 2 +(x). 



for every z G U\ \ (z 2 , Z\). 
(b) The open interval (— Zi, —z 2 ) has a neighborhood XJ 2 in the complex plane such that 



for every z G U 2 \ (—z±, —z 2 ). 

Proof. The function F = X 2+ — Ai + is purely imaginary on (z 2 , zi). Its derivative is F'(x) = 
£,2+(x) ~ = —2irip(x), and this has negative imaginary part. The Cauchy Riemann 

equations then imply that the real part of F increases as we move from the interval (z 2 , z\) 
into the upper half-plane. Thus ReA2(<z) — ReAi(z) > for z near (z 2 ,zi) in the upper 
half-plane. Similarly, ReA2(z) — ReAi(^) > for z near (z 2 ,zi) in the lower half-plane. 

By Proposition 13.21 we have ReA 3 < ReAi_ on [z 2 ,zi\. By continuity, the inequality 
continues to hold in a complex neighborhood of [z 2 , zi]. This proves part (a). The proof of 
part (b) is similar. □ 



Using the functions A^ and the constants lj, j = 1,2,3, from the previous section, we 
define 



The jump relations f)3.16j) allow us to simplify the jump matrix jx on the different parts of 
the real axis. On [z 2 , Zi], (14. 2|) reduces to 




Re\ 2 (z) < ReAi(,2) < ReA 3 (,2) 



4. First Transformation of the RH Problem 




) • (4-1) 



(4.2) 




(4.3) 



and on [—Zi, —z 2 ] to 




e n(A 2+ -Ai_) 
1 






(4.4) 
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Figure 2. The lenses with vertices —z%, —z 2 and z 2 , z\. 
On (—00, —Zi] U \—z 2 , z 2 ] U [zi, 00), (|4.2jl reduces to 

(]_ e n(A 2 +-Ai_) e n(A 3+ -Ai_)\ 
1 x G (—00, —Zi] U [—z 2 , z 2 ] U [zi, 00). (4.5) 

00 1 y 

The asymptotics of T are, because of (jl.fi |) . (jH.15|L and (j4.ll) . 

r(z) = / + 0(-) as 2 ^00. (4.6) 



Thus T solves the following RH problem: 

• T is analytic on C \ R, 

T+(x)=T-(x)j T {x), xGR, (4.7) 

• as z — >• 00, 

T(z) = J + oQV (4.8) 
Using (j4.1j) in (jl.7j) we see that the kernel K n can be expressed in terms of T as follows 

K n (x,y) = -^- -(0 e nAa +^ e nAs +W) T^^T^x) I. (4.9) 



2ni(x — y) 







5. Second Transformation of the RH Problem 

The second transformation of the RH problem is opening of lenses. Consider a lens with 
vertices z 2 ,zi, see Figure |2] The lens is contained in the neighborhood U\ of (z 2 ,zi), see 
Proposition 13.31 We have the factorization, 

/gn(Ai-A 2 ) + ^ e n(A 3+ -A!_) N 

e n(A!-A 2 )_ 

1 

(5.1 

1 W010W1 o x 

e n(Ai-A 2 )- l _ e n(A 3 -A 2 )_ | J _j g Q 1 I e ?i(Ai-A 2 ) + ^ e n(A 3 -A 2 ) + 

o o i / \ o o 1/ V o o i 
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Set 



S{z) = < 



1 

T(z) | - e "( A i(*)--Mz)) i _ e n(A 3 (z)-A 2 ( 2 )) I in the upper lens re gi OI ,. 

1 

1 

T(z) | e n (M*)-M(z)) i _ e n(\ 3 (z)~\ 2 (z)) | in the lower leng region _ 

1 



(5.2) 



Then (jOJ) and (J52J) imply that 

S+(x) = S-(x)j s (x); j s (x) 



1 N 

-l o o | , x e [z 2 ,zi] 

1, 



(5.3) 



Similarly, consider a lens with vertices —Zi, —z 2 , that is contained in U 2 (see Proposition ^. 3|) 
and set 



T(z) 



S(z) = { 



1 N 

1 ] in the upper lens region, 

_ e ?i(Ai(z)-A 3 (z)) _ e n(A 2 (z)-A 3 (z)) -y 

1 N 
T(z) I 1 ] in the lower lens region. 

e n(Ai(z)-A 3 (z))) _ e n(A2(z)-A 3 (z)) j 



(5.4) 



Then (E~7|l and (l5~l) imply that 
Set 














1 





-1 









x e [-zx, -z 2 ] 



(5.5) 



S(z) = T(z) outside of the lens regions. 
Then we have jumps on the boundary of the lenses, 

S + (z) = S-(z)j s (z), 



(5.6) 



(5.7) 



where the contours are oriented from left to right (that is, from —Z\ to — z 2 , or from z 2 to 
Zi), and where S + (S-) denotes the limiting value of S from the left (right) if we traverse 
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the contour according to its orientation. The jump matrix js in (|5.7j) has the form 

1 

j s ( z ^ = \ e n(\i(z)-\ 2 (z)) 1 e n(X 3 (z)-x 2 (z)) j on the upper boundary of the [z 2 , Zij-lens, 

1 

1 

j s (z) = ( e n ( x i( z )~ x 2( z )) i - e n(x 3 (z)-x 2 (z)) | Qn the lower boundary of the [z 2 , zi]-lens, 

1 

1 0\ 

,7s(z) = j 1 on the upper boundary of the [— zi, — 2^] -lens, 

3 n(Ai(z)-A 3 (z)) e n(X 2 {z)-X 3 (z)) j / 

1 0\ 

j's(^) = ( 1 on the lower boundary of the [— Zi, — z 2 ]-lens. 

^n(X 1 (z)-X 3 (z)) _ e n(A 2 (2)-A 3 («)) j / 

(5.8) 

On (— oo, z\] U [— z 2 , z 2 \ U [zi, oo), S has the same jump as T, so that 

S+(x) = S_(x)j s (x); j s (x) = j T (», x G (-oo, zi] U [-z 2 , z 2 \ U [zi, oo). (5.9) 

Thus, S* solves the following RH problem: 

• S is analytic on C\ (KUT), where V is the boundary of the lenses, 
• 

S + (z) = S4z) 3s (z), zGMUT, (5.10) 

• as z — ► oo, 

5(z) = / + ofiV (5.11) 



The kernel K n is expressed in terms of S as follows, see (14. 9|) and the definitions (|5.2|) and 
(|5.4j) . For x and y in (z 2 , zi) we have 

Jf n (ar,j/) = 7 r (- e nAl +^ e nA2 +^ 0) S , 7 1 (y)5 , +(z) e - nA2 +^ , (5.12) 

2m(x -y) ^ 1 

while for x and y in (— Zi, — z 2 ) we have 



K n {x,y)= e \ _ (-e" A ^fa) e" A ^)) S+\y)S+(x) I |. (5.13) 

zmyx y) \ e -nX 3+ (x) 

Since Ai+ and \ 2+ are complex conjugates on (z 2 ,zi), we can rewrite (|5.12jl for x,y £ 

(z 2 ,Zi) as 



e n{h(y)-h(x)) 

K n (x,y) = —— (_ e ™ ImA i+(f) e -namAi+(») ) ST/fj/) [ e niImAl +^ 

ZTTiyx y) 



g— ntlm Ai-|- (x)" 1 
nilm 





(5.14) 
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where h(x) = — \x 2 + ReAi + (x) as in (|1.19|) . Similarly, we have for x,y G (— z x , —z 2 ), 

n(h(y)-h(x)) / e -namAi+(x)> 

KJx,y) = - 7 (- e nilmX 1+ (y) Q g-nilm A 1+ (y) \ £-1 ,x £ + q 

2m(x-y) \ e niImA 1+ (x) 



(5.15) 



6. Model RH Problem 



As n — > oo, the jump matrix .75(2) is exponentially close to the identity matrix at every 
z outside of [—zi, —z 2 ] U [zz,Zi]. This follows from (|5.8j) and Proposition 13.31 for z on the 
boundary of the lenses, and from (|5.9|) . ()4.3|) and Proposition 13.21 for z on the real intervals 
(-00,-21), (-22,22) and (21,00). 

In this section we solve the following model RH problem, where we ignore the exponentially 
small jumps: find M : C \ ([— 21, — 2 2 ] U [22, 21]) — > C 3x3 such that 

• M is analytic on C \ ([—21, — z 2 ] U [z 2 , 21]), 

M+(x) = M-(x)j s (x), x E (-21, -2 2 ) U (22, 21), (6.1) 

• as 2 — > 00, 

M(z) = / + oQ. (6.2) 

This problem is similar to the RH problem considered in |2"U1 Section 6.1]. We also follow a 
similar method to solve it. 

We lift the model RH problem to the Riemann surface of ()2.16|) with the sheet structure 
as in Figured Consider to that end the range of the functions on the complex plane, 

n 1 =£ 1 (C\([-zi,-z 2 ]U[z 2 ,z 1 ))), 

^ 2 =6(C\[2 2 ,2i]), (6.3) 

^3=£ 3 (C\[-2i,-2 2 ]). 

Then Q x , Q 2 , f2 3 give a partition of the complex plane into three regions, see Figure EJ The 
regions fl 2 , ^3 are bounded, a 6 Q 2 , — a G f2 3 , with the symmetry conditions, 

TT 2 = n 2 , Th = n 3 , n 2 = -n 3 . (6.4) 

Denote by Tk the boundary of flfc, k = 2,3, see Figure El Then we have 

£ 1+ ([z2, z x \) = 6-([2 2 , 2!]) = r+ = r 2 n {im 2 > 0}, 
ix-{[z 2 ,z x ]) = 6+([22,2i]) = r 2 - = r 2 n {im2 < 0}, 
£i+([-*i, -22]) = e 3 -([-2i, -2 2 ]) = r+ = r 3 n {im2 > 0} 
6-([-2i, -2 2 ]) = -2 2 ]) = r 3 = r 3 n {im 2 < 0}. 

We are looking for a solution M in the following form: 

/Mi(6(2)) Mi(6(2)) Mi(6(z))\ 
M(z) = M 2 (6(2)) M 2 (6(2)) M 2 fe(2)) , (6.6) 
\M 3 (a(2)) M 3 (6(2)) M 3 (6(2))/ 
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a 1 




f a 3 ) 


^ — ^\ 

f Q 2 ) 


-q[ -a J-p 


p[ a U " 







Figure 3. Partition of the complex £-plane. 

where Mi(£), M 2 (£), M 3 (£) are three analytic functions on C \ (r\ U T 2 ). To satisfy jump 
condition (|6.1j) we need the following relations for k = 1,2, 3: 



ii4 + (£) = M fc _(£), £er 2 -ur 3 , 
M fc+ (0 = -M fc _(0, ^r+ur+. 

Since £i(oo) = oo, £2(00) = a, £,3(00) = —a, then to satisfy (|fi.2j) we demand 

Mi (00) = 1, Mi (a) = 0, Mi(-a) = 0; 
M 2 (oo) = 0, M 2 (a) = 1, M 2 (-a) = 0; 
M 3 (oo) = 0, M 3 (a) = 0, M 3 (-a) = 1. 

Equations ()6.7j) - (j6.8j) have the following solution: 

e-a 2 „ , , e±« 



(6.7) 



(6.8) 



Mi(0 



M 2 , 3 (0 = C 2> ; 



(6.9) 



V(e-P 2 )(e-? 2 ) 'VF?)F? 

with cuts at r 2 , r 3 . The constants c 23 are determined by the equations M 23 (±a) = 1. By 

(Q, 

hence 



" P 2 )(£ 2 - g 2 ) = £ 4 - (1 + 2a 2 )^ 2 + (a 2 - l)a 2 , 

2a 



M 2 a 



c 2 - 



By taking into account the cuts of M 2 (£) we obtain that 

M 2 (a) = c^v^, 

hence 

% 



(6.10) 
(6.11) 

(6.12) 
(6.13) 
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Similarly, 



hence c 3 is the same as c 2 , 



M 3 (-a) = c 3 



c 3 



-2a 



c 3 « 



V2, 



(6.14) 



(6.15) 



Thus, the solution to the model RH problem is given as 



M(z) 



• 6jg)±g 



V«3(*)-P 2 )«30O-9 2 ) 



^(^CO-p 2 )^*)-* 2 ) 

\ V^W-P 2 )^)-? 2 ) ^2(^f( 2 )-p2)(^W- ? 2) V 2 «3( 2 )-P 2 )«lW-9 2 )/ 



^w^iw V 2 (^)-p 2 )(e 2 (z)-q 2 ) 

6(g)-g „• &2(z)-a 



— Z 



—2 



(6.16) 



with cuts on [z 2 , Zi] and [— zi, —z 2 }. 

The model solution M(z) will be used to construct a parametrix for the RH problem for 
S outside of a small neighborhood of the edge points. Namely, we will fix some r > and 
consider the disks of radius r around the edge points. At the edge points M(z) is not analytic 
and in a neighborhood of the edge points the parametrix is constructed differently. 

7. Parametrix at Edge Points 

We consider small disks D(±Zj,r) with radius r > and centered at the edge points, and 
look for a local parametrix P defined on the union of the four disks such that 

• P is analytic on D(±Zj, r)\(lU T), 



P+(z) = P-{z)j s {z), ^(MUf)n D{± Zj ,r), 



(7.1) 



as n — > oo, 



P(z) =[1 + 



n 



M(z) uniformly for z G dD(±Zj,r). 



(7.2) 



We consider here the edge point Z\ in detail. We note that by (|3.11|) and (|3.14j) we have 
as z — > Zi, 

2pi 



so that 



X 1 (z) = q(z - Zl ) + -f(z - Zl f 2 + 0(z - Zl f 
\ 2 {z) = q(z - z x ) - ^(z - Zl f' 2 + 0(z - z x f 

X 1 (z) - X 2 (z) = ^(z - Zl f/ 2 + 0(z - Zl f' 2 



as z — > z\. Then it follows that 



P{z) 



-(X^z) - X 2 (z)) 



2/3 



(7.3) 



(7.4) 



(7.5) 
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is analytic at z\, real- valued on the real axis near z\ and /3'(zi) = p^ 3 > 0. So (3 is a 
conformal map from D(zi,r) to a convex neighborhood of the origin, if r is sufficiently small 
(which we assume to be the case). We take T near z\ such that 

(3(T n D{z x ,r)) C {z | arg(z) = ±2tt/3}. 

Then T and K. divide the disk D(z\,r) into four regions numbered I, II, III, and IV, such 
that < arg (3{z) < 2tt/3, 2n/3 < axg/3(z) < tc, -ix < arg (3{z) < -2n/3, and -2tt/3 < 
axg/3(z) < for z in regions I, II, III, and IV, respectively. 

Recall that the jumps j$ near z\ are given by (|5.3|) . ()5.8|) . and ()4.3|) : 



J 5 



Js 



Js 




on zj - r,zx J 



1 

e n(Ai-A 2 ) ]_ e n(A 3 -A 2 ) on U pp er boundary of the lens in D[z\ , r) 







e n(Ai-A 2 ) y 






(7.6) 



on the lower boundary of the lens in D(z\,r) 



on (zi, zi+r\. 



We write 



P 






1 



in regions I and IV 
in regions II and III. 



(7.7) 



Then the jumps for P are P + = P-jp where 

jp = | -i u u | on [zi - r,zi) 

Jp 




1 > 

e n(A!-A 2 ) l o on the upper side of the lens in D(z±, r) 







1 



Jp 



Jp 



1 > 

g n(Ai-A 2 ) i o I on the lower side of the lens in D(zi, r) 

1 / 

'\ e n(A 2 -Ai) gN 

1 I en ( :,. ■■: /-|. 

,0 1, 



{IX 
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We still have the matching condition 



P( z ) = ( I + O ( - ) J M(z) uniformly for z e dD(z u r) 



(7.9) 



since Re A3 < ReA2 on D(z%,r), which follows from Proposition 13.21 

The RH problem for P is essentially a 2 x 2 problem, since the jumps (J7.8j) are non-trivial 
only in the upper 2x2 block. A solution can be constructed in a standard way out of Airy 
functions. The Airy function Ai(z) solves the equation y" = zy and for any e > 0, in the 
sector 7r + e < aigz < it — e, it has the asymptotics as z — » 00, 



Ai(2) = v^ e " , " > ' 2 ( 1+0(2 " /2) > 



(7.10) 



The functions Ai(uz), Ai(tu 2 z), where lu = e 2 ™, also solve the equation y" = zy, and we 
have the linear relation, 



Write 



M(z) + uok\(ujz) + uo 2 Ai(uj 2 z) = 0. 



y (z) = Ai(z), yi(z) = uAi(uz), y 2 {z) = u 2 Ai(u 2 z), 



(7.11) 
(7.12) 



and we use these functions to define 



$(z) = < 



'y (z) -y 2 (z) N 
y' (z) -y' 2 (z) 
1, 

- yi {z) -y 2 (z) N 
-y{(z) -y' 2 (z) 
I 

'-y 2 (z) y x {z) N 

01 



for < argz < 2n/3, 

for 27r/3 < argz < 7r, 
for — 7r < argz < — 27r/3, 



(7.13) 



y' (z) y[(z) 0|, for -27r/3 < argz < 0. 







1 



Then 



P(z) = E n (z)$(n 2/3 /5(z))diag (e^ n{Xl e -^(AiW-A a (*)) j j 



(7.14) 



where _E n is an analytic prefactor that takes care of the matching condition (|7.9|) . Explicitly, 
is given by 



1 -1 N 
£ n = v/vrM I -z -i 



■n 1 / 6 ^ X 
o 
1 



(7.15) 



1 / 

A similar construction works for a parametrix P around the other edge points. 
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ID 



Figure 4. The contour T R for R. 

8. Third transformation 

In the third and final transformation we put 

R{z) = Siz)M(z)- 1 for z outside the disks D(±zj, r), j = 1, 2 

R(z) = S(z)P(z)~ l for z inside the disks. 

Then R is analytic on C \ Tr, where Tr consists of the four circles dD(±Zj,r), j = 1,2, 
the parts of T outside the four disks, and the real intervals (— oo, —z\ — r), (— z 2 +r, Zi — r), 
(zi + r, oo), see Figure |U There are jump relations 

R + = R^j R (8.2) 

where 

jr = MP^ 1 on the circles, oriented counterclockwise 
]r = MjgM^ 1 on the remaining parts of Tr. 

From ()7.2|) it follows that Jr — I + 0(l/n) uniformly on the circles, and from ()5.8|) . ()5.9|) 
()4.3|) and Propositions 13 .21 and 13. 31 it follows that ]r = I + 0(e~ cn ) for some c > as n — > oo 
uniformly on the remaining parts of Tr. So we can conclude 



■3) 



n 



as n — > oo, uniformly on Tr. 



■4) 



As ^ — > oo, we have 

R(z)=I + 0(l/z). (8.5) 
From ()8.2|) . (|8.4j) . ()8. 5|) and the fact that we can deform the contours in any desired 
direction, it follows that 

1 



R( z ) =1 + 



n(\z\ 



as n 



oo. 



(8.6) 



uniformly for z E C \ T R , see O CS1 d d] • 

By Cauchy's theorem, we then also have 

R'(z) = O 



n(\z\ + 1) 
and thus 

R-\y)R(x) =1 + R-\y) (R(x) - R(y)) 
which is the form we will use below. 



I + O 



x-y 



n 



17) 
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9.1. Proof of Theorem 11.11 Consider x G (z 2 ,zi). We may assume that the circles 
around the edge points are such that x is outside of the four disks. Then ()8.1|) shows 
that S(x) = R(x)M(x) and it follows easily from ()8.7|) and the fact that M + is real analytic 
in a neighorhood of x that 



S + 1 (y)S+(x) = I + 0(x-y) as y 
uniformly in n. Thus by (|5.14|) we have that 

n(h(y)-h(x)) 



X 



(9.1) 



K n (x,y) 



2iri(x — y) 

e n(h{y)-h(x)) 
n(h(y)-h(x)) 



g— nilm Ai+ (ie) n 

f e nilmXi+(y) g-nilm Xi+(y) q) (J + O (x — y)) I e " iIm -^i+M 



_ e m(Im Ai + (j/)-Im \\+{x)) _|_ g — ni(Im Ai+(2/)-Im \\ + (x)) 



2ni(x — y) 
sin(nlm (A i+ (x) - \ 1+ (y))) 





0(1) 



vr(x - y) 



+ 1 



and the 0(1) holds uniformly in n. Letting y — > x and noting that by (jH.14j) and ()H.9)1 



(9.2) 



— ImX 1+ (y) = lm£ 1+ (y) = np(y) 



we obtain by l'Hopital's rule, 



K n (x,x) = np(x) + 0(1), 



(9.3) 



(9.4) 



which proves Theorem 11.11 if x G {z2,Z\). The proof for x G (—Zi, —z 2 ) is similar, and also 
follows because of symmetry. 

For x G (— oo, — Zi) U (— z 2 ,2; 2 ) U (zi,oo), we have that K n (x,x) decreases exponentially 
fast. For example, for x > z±, we have that 



K n {x,x) = O (e- n{Xlix) ' X2{x)) ) 



as n — > oo. 



(9.5) 



This follows from ()4.9|) and the observation that that T + 1 (y)T + (x) = I + 0(x — y) as y — >• x 
if x > ^i- It is clear that ()9.5|) implies 



1 

n— >oo 



lim —K n (x, x) = 0. 



(9.6) 



We also have (|9.6|) if x is one of the edge points. In fact, for an edge point x it can be 
shown as in the proof of Theorem 11.31 that 



n 



-K n (x,x) = O 



n 



1/3 



asn-> oo. 



This completes the proof of Theorem 11.11 



(9.7) 
□ 
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9.2. Proof of Theorem 11.21 We give the proof for xq G (z2,Z\), the proof for xq G 
{—Zi,—Z2) being similar. We let 

u v 

x = x ^ — r, y = x H -, — r. (9.8) 

np(x ) np{x ) 

Then we have ()9.2j) . and so by the definition p. 18)1 of K n , 

1 -k n (x, y) = Sin(nIm (A1 ; (X) ~ Xl+iy))) +o(-). (9.9) 



np(xo) ' ti(u — v) \n j 

Because of (I9.3J1 we have by the mean value theorem, 

Im (Ai+(x) - A 1+ (y)) = (x - y)Trp(t) (9.10) 
for some t between x and y. Using ()9.8j) we get t = xq + 0{l/n) and 

nlm {X 1+ {x) - Xi+(y)) = n(u - v)-^ = n(u -v)(l + o(^j. (9.11) 

Inserting (|9.11|) into (j9.9j) . we obtain 

1 a . . sin7r(u — i>) „ / 1 \ . 

K n (x, y) = — ^ ^ + O - (9.12) 



np(x ) ' 7r(u — u) yn 

which proves Theorem 11.21 □ 

9.3. Proof of Theorem II. 3L We only give the proof of (jl.22j) . since the proof of (jl.23j) is 
similar. We take p\ as in ()3.11|) and we recall that (3'{z\) = p 2 / 3 . 
Take u, v G K and let 

u v 

x = z i + 7 w73> y = z i + 7 ^77- ( 9 - 13 ) 

(pin) 2 '- 1 (pin) 2 ' d 

Assume u, v < so that we can use formula (|5.14j) for K n (x, y). Then we have that x belongs 
to D(zi,r), for n large enough, so that by (|8.1jh ()7.7|) and ()7.14|) 

= i?(x)P + (x) = R(x)P+{x) 

= R(x)E n (x)$ + (n 2/3 p(x))diag (e^-^O+O*), e -|n(Ai-A a )+(*) ) ^ ( 9 . 14 ) 
= J R(x) J E„(x)$ + (n 2/3 /3(a;))diag ( e «i^i+W )e -" iIm ^W ) x ) 
and similarly for S + (y). Then we get from (|5.14|) and (jl,18|) 



(pin) 2 / 3 n ' 2m(u — f ) 

'1\ (9.15) 
x i?(x) J E„,(x)$ + (n 2/3 /5(x)) I 1 



.0, 



Since p 2//3 = (3'(zi), we have as n — > oo, 



n 2 / 3 /3(x) = n 2 / 3 /3 ( + — -p j - it (9.16) 
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which implies that $ + (n 2//3 /3(x)) — ► We use the second formula of (j7.13j) to evaluate 

(since u < 0), and it follows that 

-yl(tt)-^(tt) I I .'/,(") 1 • ( <: ^!') 

Similarly 

lim (-1 1 0) ®^{n 2, *p(y)) = (-1 1 0) §+\v) 



(9.18) 

= -2iri (y' 2 (v) + y[(v) -y 2 (v) - y x {v) 0) 
= -2m(-y' (v) y (v) 0) . 

The factor —2iri comes from the inverse of $+(?;), since det $ = (— 2-ni)~ l by Wronskian 
relations. 

Next, we recall that R- 1 (y)R(x) = 1 + (2=«), so that by (f9~T3j) 

J R- 1 (y)i?(x)=/ + of z ^V (9.19) 





(yo{u)\ 


)-l 


y' (u) 1 








2/2 («) 




(v) 









(9.20) 



The explicit form ()7.15|) for E n readily gives 

E n (x) = 0(nV»), ^(y) = 0(n l /«), i?" 1 (?/)£„(*) = / + O 

Combining ()9.19|) and (|9.20|) . we have 

lim E-\y)R-\y)R(x)E n (x) = I. (9.21) 

n— >oo 

Inserting (jHUZJ), (f9~T%|) . and (JOT]) into fEHKl) . we obtain 



1 



lim -^ K n&y) = tt-( v x (- 27ri ) (-y'o( v ) yo( v ) °) y' {u) 

n^oo cn ziri{u — v) V 7 (9 22) 

_ yo{u)y' {v) - y' (u)y (v) 



u — V 

Since y — Ai, we have now completed the proof of (|1.22|) in case u, v < 0. 

For the remaining cases where u > and/or t> > 0, we have to realize that we have not 
specified the rescaled kernel K n (x,y) for x and/or y outside of [—Zi, —z 2 \ U [22, zi], since in 
(I1.19|) h is only defined there. We define 

1 1 

h(x) = --x 2 + - (Ai(x) + \ 2 (x)) , xe(z!,oo). (9.23) 

We will assume in the rest of the proof that u > and v > 0. The case where u and v 
have opposite signs follows in a similar way: then we have to combine the calculations given 
below with the ones given above. 
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So let u, v > and let x and y be as in ()9.13j) . For the kernel K n we start from the 
expression ()4.9|) in terms of T. Since u > 0, we have x, y > Z\, and so we have by ()5.6|) . 
(1771) and (EH, 

/l 

T+(a;) = S+(x) = R(x)P + (x) = R(x)P+(x) 1 e "(M*)-W) 

\0 

= J R(x)E n (x)$ + (n 2/3 /3(x)) (9.24) 

/ e |n(Ai(*)-Aa(a!)) Q 0\ /l 

X I e -in(AiOr)-A 2 (V)) Q J I Q 1 e «(A 3 W-A 2 W) 

V 1/ \o 1 

Then 

r+(ar) I j =e-^ n{Xl(x) - X2{x)) R(x)E n (x)^ + (n 2/3 f3(x)) I 1 . (9.25) 

As before, we have $ + (n 2//3 /3(x)) — ► ( &+(u) as n — ► oo. Now we use the first formula of (|7.13j) 
to evaluate $+(«) so that 

A 



T 

lim $+(n 2/3 /5(x) [ ] = $+(«) 
,0. 



2/o W 

1 = 1 I • (9-26) 

V / V 



We have as in 

/l \ / e -3"( A i(if)--MiO) o 0\ 

T" 1 ^) =01 -e< X ^~ X ^ e in(A 1 fe)-A 2 (j / )) q 

\o o i /V o o 1/ ( 9 - 27 ) 

x $; 1 (n 2 / 3 /3(y))E- 1 (y) J R- 1 ( 2/ ), 

so that 

(0 e nX ^ e nX3 ^)T- l (y) = e 1 2 n ^- x ^ (0 1 0) ^- 1 (n 2/3 (3(y))E~ 1 (y)R- 1 (y). (9.28) 
We have 

lim (0 1 0)^l 1 (n 2/3 (3(y)) = (0 1 0) $7 1 (m) 

n-+oo 2g^) 

= {-2m) (-y' (u) y (u) 0) 

and as before we have (|9.21|) . 

Inserting (ET231) and (ET2"K1) into and using the limits (j9~2Tj) . (ET2T)j) and (ET2HJ) . we 

arrive at ()1.22|) in the case u, v > 0. 

This completes the proof of Theorem 11.31 □ 

10. Large n asymptotics of the multiple Hermite polynomials 

As noted in Section |21 see also [0], we have that the (1, 1) entry of the solution Y of the 
RH problem (jl.4j) - (jl.5|) is a monic polynomial P n of degree n satisfying 

P n (x)x k Wj(x)dx = 0, k = 0, 1, . . . , rij — 1, j = 1, 2. 
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For the case Wi(x) = e - n (h x2 ~ ax \ w 2 (x) = e - n i\ x2 + ax ) ( this polynomial is called a multiple 
Hermite polynomial The asymptotic analysis of the RH problem done in Sections 

EHH1 also yields the strong asymptotics of the multiple Hermite polynomials (as n — > oo with 
n even and n\ = n 2 ) in every part of the complex plane. We describe these asymptotics 
here. Recall that P n is the average characteristic polynomial of the random matrix ensemble 
ffTTjl . see (Q|) . 
We will partition the complex plane into 3 regions: 

(1) Outside of the lenses and of the disks D(±Zj,r), j = 1, 2. 

(2) Inside of the lenses but outside of the disks. 

(3) Inside of the disks. 

We will derive the large n asymptotics of the multiple Hermite polynomials in these 3 regions. 
(1) Region outside of the lenses and of the disks. In this region, we have by (|5.6J1 and 

T(z) = R(z)M(z), (10.1) 

hence by f|4.1|) 

diag(e- n \ e~ nh , e- nl3 )Y{z)di&g{e-^ z2 , e~ naz , e naz ) 



(10 2) 

= J R(z)M(z)diag(e" nAl(2) , t~ nX ^ z \ e""^). 
By restricting this matrix equation to the element (1, 1) we obtain that 

3 

P n (z)e-^ z2 = e'^J^^Wjiiz), (10.3) 

where 

A(*) = AiCz)-Zi = J £i(s)ds, (10.4) 

and as z — > oo, 

z 2 

\{z) = — -\nz + 0{z- 2 ). (10.5) 
In the sum over j in ()10.3j) the term j = 1 dominates and we obtain because of (joMfij) that 

where for the square root we use the principal branch (the one that is positive for z > z\), 
with two cuts, [—2i, — zj\ and [z%, Z\]. 

(2) Region inside of the lenses but outside of the disks. In this region, we get from 
Q , (fOjl and dHUJ), 

T(z) = R(z)M(z)L(z)- 1 , (10.7) 
where L(z) is the matrix on the right in (|5.2jl and (|5.4jl . Hence by (j4.1j) 



n _2 



diag(e-^,e-" (2 ,e-" (3 )r(z)diag( e -^ , e -" az ,e" a2 ) 

= J R(z)M(z)L(z)- 1 diag(e- nAl(2) ,e- nA2(2) ,e- nA3(z) ). ' 
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Consider z the upper lens region on [z 2 , zi\. Then 



1 

Liz) = I — e n ( Al (*)~ A2 ( z )) 1 — e «( A 3(2)-A 2 (^)) 
1 



(10.9) 



hence 



L(z 



,-i 



1 

e n(Ai(z)— A2O2)) ]_ e n(A 3 (2)-A 2 (z)) 
1 



(10.10) 



and the first column of the matrix M(z)L(z) 1 diag(e n M z ) ; e nA 2( 2 ) )e iA 3 («)^ j g 

/ M 1 (^ 1 (2))e- nAl ^ + M 1 (^ 2 (2))e~ nA2 (^\ 

M 2 (ei(^))e- nAl ( 2 ) + M 2 (&(z))e- nX2 ^ , (10.11) 
v M 3 (ei(^))e-" Al ( 2 ) + M 3 (&(z))e- nX ^ J 

see (16. 6|) . By restricting equation ()10.8|) to the (1, 1) entry, and using ()6.16j) and ()8.6|) . we 
obtain that in the upper lens region on [z 2 , Zi] 



g(z) - a 2 



el(^) - a 2 



+ 



+ 



oil 



-nXi (z)+nl\ 



where 



M 2 ) = / rfs > fc = 1,2. 



(10.12) 



-nA 2 (z)+nZi 



(10.13) 



In the same way we obtain that in the lower lens region on [Z2, z%], 



PJz)e-^ 2 



+ o(- 

n 



-n\i (z)+nli 



+ o(- 

n 



(10.14) 



-nA 2 (z)+n/i 



For z = x real, x 6 \zi + r, zi — r] , both (|10.12|) and ()10.14j) can be rewritten in the form 

1 



where 



and 



A(x) cos[nIm Xi + (x) — <p(x)] + O 

i 2 l+ (x)-a 2 



n 



-n Re Ai+ (x)+nl\ 



A(x) = 2 



<f(x) = arj 



(10.15) 



(10.16) 



(10.17) 
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By using equation ()H.9j) . we can also rewrite (|l().15jl in terms of the eigenvalue density 
function p(x), 

jl(a;)cos nir £ p(s) ds - <p{x) + O J e -«R**i+(*)+"*i. (10.18) 

Equation (|l().18jl clearly displays the oscillating behavior of P n on the interval [z 2 + r, z\ — r\. 
It also shows that the zeros of P n (x) are asymptotically distributed like p(x)dx, the limiting 
probability distribution of eigenvalues. Similar formulae can be derived on the interval 
[—Zi + r, —z 2 - r]. 

(3) Region inside of the disks. Consider the disk D(zi,r). In the regions I and IV, we 
have by JZZZD, (EJ) and flESD 

/l 

T(z) = R{z)P{z) = (I + 0{n- 1 )) P{z) 1 e n ( A3 ^- A2 ^) | , (10.19) 

\0 1 

hence by (@~TJ), (T7TTD . and (T7TTHJ) 

diag(e~"\ e^ 2 , e- nZ3 )r(z)diag(e-t 22 , e^, e na ") 

/ T2 1 /6 / 5( 2 )l/4 _ n -l/6 /9 ( 2 )-l/4 Q N 

= (/ + 0(n- 1 )) v ^FM(^) -in 1 / 6 /?^) 1 / 4 -m" 1 ^^)-!/ 4 

\ 1/ (10.20) 

/l N 

x $(n 2/3 /?(z))diag(e- nQ(2) ,e- na(z) ,e- nA3(2) ) 1 1 

\0 1. 

where 

a(z) = Xl{z) + 2 HZ) . (10.21) 

By restricting equation ()10.20|) to the (1,1) entry, and using the first expression of (|7.13j) 
(in region I) or the fourth expression of (J7.13j) (in region IV) to evaluate &(n 2 / 3 fl(z)), and 
()6.16|) to evaluate M(z), we obtain that 



P n {z)e-% z = v^F [n 1/6 B(z)Ai(n 2/3 P(z))(l + 0(rr 1 )) 

+n- 1/6 C(z)Ai'(n 2/3 (3(z))(l + 0(n~ 1 ))] e~ na{z)+nh 



10.22) 



where 



D( ) B(z) 1/4 [ ^ = ^ - i ^(^)-« 2 | , () 

' Vfe 2 W-p 2 )(^)-g 2 ) V^W-P 2 )^ 2 ^)-? 2 ) 



and 



C(z) = /?(^)- 1/4 ( - - ^)~° ! j - ^ 2 ^)~ a2 ] . (10 .24) 
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The same asymptotics, p().22jl . holds in regions II and III as well. Thus, (110. 22j) holds in 
the full disk D(zi,r). It may be verified that the functions B(z) and C(z) are analytic in 
D( Zl ,r). 

This approach allows one to derive a formula similar to (|10.22|) in all the other disks 



D(±zj, r) as well. 



Appendix A. Recurrence equations for multiple Hermite polynomials 
From orthogonality equation (|2.2|) . we obtain that as z — > oo, 



1 

2vri 



P ni , n2 (u)w k {u) 



du 



u — z 



2m 



P nit n 2 (u)w k (u) ( - + + 



z Z^ 



1 / h( k ) (fc) 



2m \ z nk+1 z nk+2 



du 



fc=l,2, 



(a.i; 



where for k — 1, 2, fci£ )na is defined in (|2.Hjl and 



q {k) 



This implies that 



*ni,n a (*) 



/ + 



P niin2 (x)x nfc+1 w fc (a;)da;. 



diag ^z n e 



n -\Nz 2 —l —ni„-Naz „—\„—n2„Naz\ 



Cn Z 



i u 2 



(A.2) 



(A.3) 



where 



ni,n 2 



V 





ft' 2 ' 


"ni — l.no 
L 'n 1 -l,n 2 


h (2) 

"ni ,719 — 1 
"n 1 -l, 112 


u n 1: n 2 -l 


(2) 




/i (2) 



/ + 



(1) 



+ 



/ 2 

z 



V 



o 



2P, + v]> (1) , A - P,«r( 1 ) + 
"IT ^ni+l,ri2 1 1 ni,n 2 ~ 



-1 ft ni,n.2 





/l (2) 

' t, ni+l,n2— 1 
lb n 1 ,n2 — l 



\ 





,( 2 ) 



* 711,712 



(A.4) 



and P nitn2 (z) = z n + p ni ,n 2 z n 1 H • Set 

Then by (|2.9|) . C/ niin2+ (x) = f/„, 1>n2 _(x) (i.e., no jump on the real line) and as z — > oo, 



(A.5) 



+ 



+ 1 ,ri2 —1 



P 3 + 



where 



Pi = diag(l,0,0), P 2 = diag(0,l,0), P 3 = diag(0, 0, 1) 



(A.6) 
(A.7) 
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Since U ni}Jl2 (z) is analytic on the complex plane, equation (|A.6|) implies, by the Liouville 
theorem, that 

,(2) 



hi 

'3 

711 ,712 — 1 



TJ ( z \ _ ~D , ^(1) D _ _ (_ "l + l,"2-lp f 

<^ni,ri2l^J ' J 1T ^ni+l,n2 J 1 1 1^711,712 ' , (2) 1 

ft, ' 

% "711,712 Cni,ri2 ^711,712 

1 

1 e 

1 w °ni,n2 



(Ai 



where 

h (1) h (2) h {2) 

'mun2 / n J 'tni,n 2 , „ ' t, ni+l,n 2 -l / ^ / * q\ 

c ni,n2 /-n T 1 "> "ni,n2 / 2 ) ' ' e 7ii,ri2 (2) ' v.-™""/ 

711 — 1,712 711,712—1 711,712 — 1 

Thus, we obtain the matrix recurrence equation, 

^711+1,712 (^0 ^711,712 (^0^711,712 

(2;). (A.10) 

By restricting it to the element (1, 1) we obtain that 

-^711+1,712 ^711 ,712 ) 711 ,712 (^0 '-'TO ,712 -^711 — 1,712 (^0 ^711,712-^711,712— l(^0 ) (A. 1.1. J 

and by restricting it to the element (3, 1) we obtain that 

1,712 ~ 1 (^0 -^711,712 (^0 "I" ^7ll ,712 -f Tli ,"2 — 1 (^0 ' (A. 1.2) 

Similar to (|A.10|I . we have another recurrence equation, 

*m,n2+l(^) = ^ni.na^)*^^^), (A.13) 

where 

^711,712 < -' Jlli?1 2 "TO,712\ 

1 e m>n2 , (A. 14) 

1 ) 

and 

tfl-l n 2 +l 

= jjy 2 7^ 0- (A. 15) 

^711-1,712 

By restricting (|A.13j) to the elements (1, 1) and (2, 1), we obtain the equations, 

-Prii,n2+l(^0 ^71 1 ,712 ) -^711,712 (^0 ^m,n2 -P711— 1,712 (^J ^711,712-^711,712— l(^) > (A. 16) 

and 

-P7II — 1,712+1 V^) -^"7il,n2(^0 "I" ^711 ,712 Pn\ — 1,712 (^0 ' (A. 1.7) 
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Appendix B. Differential equations for multiple Hermite polynomials 
Set 



An 1 ,n 2 {z)=xK 1 ,n 2 {z)*n 1 ,n 2 {z) (B.l) 

It follows from ()2.9|) . that A ni>n2 (z) has no jump on the real axis, so that it is analytic on 
the complex plane. By differentiating (jA.3|) we obtain that as z — > oo, 



A ni ,n,2 (^) 



/ + 



+ 



o cr 

-a | | / 

a, 



+ 



z 



(B.2) 



Since y4 niin2 (z) is analytic, we obtain that 

+ • • 



^4ni,n 2 (^) 



(<* 






2! 









—a 





!■ 










(T 
0|(/ 

.0 0, 



+ 



pol 



(B.3) 



where [/(z)] po i means the polynomial part of f(z) at infinity. From (|B.1J) we get the differ- 
ential equation, 

K u nM = NA ni>m {z)V nun2 {z). (B.4) 

and (IB .3D reduces to 



"2 £-ni,7i2 dni,n2 

-I -a 
-10 a 



(B.5) 



Appendix C. Proof of Proposition 12. f I 

Equations (|A.10|I . ()A.13|) . (jB.4|) form a Lax pair for multiple Hermite polynomials. Their 
compatibility conditions are 



1 .-. 



N 



U nil n 2 ( Z ) — A ni ^ n2+ i(z)U ni ^ n2 (z) — U nim (z)A ni>n2 {z). 



This gives the equations, 



(CI) 



^ni,n 2 ^ni+l,fi2 '-711,712 "I - > ( ^ni+l,n 2 ^ni,ri2) 6 niin2 2(2, 



"711,712 O) C reijn2 -I-1 ^7ii,Ti2> 

Since co, n2 = ^ni,o = 0, we obtain that 



+ 1 



°ni,n 2 j "711,712 



dni,U2 yy) ^711,712 2fl- 



n2 
AT 



(C.2) 



(C.3) 



This proves the first equation in (|2.10j) and equation (|2.11jl . Similarly we obtain that e„ 1)n2 = 
2a and this proves the second equation in ()2.10|) . Proposition 12. II is proved. □ 
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